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Abstract 

We present a new expression of the scattering amplitude, valid for spher- 
ical absorbing objects, which leads to an improved version of the eikonal 
method outside the diffraction region. Limitations of this method are dis- 
cussed and numerical results are presented and compared successfully with 
the Mie theory. 
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1 Introduction 



The eikonal picture, or eikonal model, has been first developped by Perrin 
and Chiappetta [Q to describe optical scattering phenomenons by scatter- 
ers having an axis of symmetry parallel to the incident wave- vector. This 
approximation is a variant of the High-Energy Approximation (HEA) in- 
troduced by Glauber for scattering of particles by a potential of infinite 
range. Even though this latter theory is valid only for small scattering angle, 
the eikonal picture was proven to be valid in the whole angular domain. Nu- 
merical results [Q-Q for the sphere shown rather good agreement with the 
Mie theory, particularly when internal multiple refiexions are negligible. This 
useful approximation provides a simple mathematical scheme and has been 
used particularly in astronomy, to describe the light scattering by irregular 
particles [^. 

An other version of the eikonal method, called the generalized eikonal 
approximation including the second Born correction of the scattering 
amplitude has been developped for the sphere, leading to good agreement 
with the Mie theory. Recently, a vectorial version of the eikonal method 
valid for objects of approximately spherical shape has been derived 0. The 
eikonal approximation has also been extended to the case of objects whose 
axis of symmetry is not parallel to the direction of the incident wave-vector 

i- 

The eikonal picture applied to scattering by large spheres gives only a 
qualitative description of the backward intensity. In the present work, we 
show that simple considerations suggest a modification of the eikonal model 
amplitude for large absorbing spheres, leading to a better quantitative de- 
scription of the scattering phenomena. 

This paper is organized as follows : in section 2, we start by recalling the 
eikonal model and we give the scattering amplitude for an absorbing sphere. 
In section 3, we consider this scattering amplitude in the forward and back- 
ward limits. Comparisons with the diffraction and the optical backscattering 
amplitudes suggest an empirical transformation of the scattering amplitude 
which is tested against exact solutions in section 4. The conclusion is given 
in section 5. 
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2 The eikonal scattering amplitude for an ab- 
sorbing sphere. 

The eikonal method leads to an approximate formulation of the solution of 
the scalar Helmholtz equation 

{a + m2A;2}^(r) = 0, (1) 

where m is the refractive index of the scatterer, k is the wave-vector of 
the incident plane wave, and \E'(r) is the electric field perpendicular to the 
scattering plane For objects having an axis of symmetry parallel to k the 
solution \E'(r) can be written for large r : 

^(r)=e^'^- + /(^)— + 0(i). (2) 
The scattering amplitude is given in cylindrical coordinates, by : 

f(9) = / dbbJoikb sin e)G{b) + 0(1), (3) 

2 JO 

where G{b) is 

dzUib,z)e-''''exp{i- dz'Uib, z')}, (4) 

z{b) is the boundary of the scatterer for a given value of the impact parameter 
b, and U{b, z) = 1 — rrP, q = 2A;sin^ |. The incident wave-vector k has been 
taken parallel to the z axis (see Fig. 1). The function G{b) can be computed 
exactly when the refractive index is a constant. For a convex particle whose 
center of symmetry is located at the origin : 

G{b) = '^^{e){e-"^'^'^+'^^''^ - (5) 

K 

kU 

with 7(6*) = ; — and the eikonal function x(6) reads : 

2q — kU 

k r^W 

Xib) = - dz'U{b,z'). (6) 

A J— 00 
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For a sphere of radius a, x(6) = kUz{h). Then the scattering amphtude has 
the simple expression : 

f{e) ^ -ik-f{e) £ dbbJo{kb sin 9) [6'''^'^^^+''''^^'' - e''^'^^^}, (7) 



with z{b) = — b"^- 

When the refractive index has a non neghgible imaginary part, i.e. when 
the condition \Im{'m)\ka > 1 is satisfied, the real part of ix{b) takes a high 
negative value which leads to the approximate expression of f{0), namely : 

f{e)^ik-f{e) r dbb Jo{kb sine) e^^'^^l (8) 
Jo 

3 Modification of the eikonal model 

Assuming that the sphere is sufficiently absorbing means that physically, 
the diffractive and the first order reflective parts of the scattering amplitude 
give a reasonable description of the scattering phenomenon. The present 
expression (|^) of the scattering amplitude leads only to qualitative results 
in the backward direction and needs to be modified. Let us first consider 
the behaviour of f{9) in when 9 is small ( i.e. when the diffractive part 
dominates). Since, e*'^^*^^-' ~ 1, we obtain : 

fi9)^--J,{kasm9), (9) 
sm U 

showing that we recover the exact expression of the diffraction amplitude for 
a sphere. On the other hand, in the neighbourhood of the backscattering 
angle ( i.e. |7r — 6*1 << 1), the use of the method of stationary phase leads 
to : 

f(9) ^ 

2 sin' I " ' 2 



/(^) ^ 77— ^exp(2^fcasin2-)7(0). (10) 



Since sin 2 is close to 1, we approximate (|TOD as 



/(^)«|exp(2^fcasin^)7(^). (11) 



Whereas the geometrical optics gives : 



CL 9 

fopt ~ -exp (2iA;asin -)r^{9), (12) 
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where r±{6) ^ is the Fresnel reflection coefficient for the perpendicular com- 
ponent of the electric field, which was proven in [Q] and to be identical to 
the eikonal amplitude. We observe that the expression ( ]TT| ) differs from (|I^) 

by a factor — -— which has a modulus different from 1 for refractive index 
r±i9) 

near 1, leading for large angle to a systematic error in the evaluation of the 
scattering intensity. However, since r±{0) = 7(0), a simple modification of 
the eikonal model which takes into account the forward and the backward 
asymptotic behaviors (||) and ([T2|) , would consist in replacing 7(6') by r^{9), 
leading to : 

fmodiO) ~ ikr^ie) r dbbJoikb sin e)e'^''^'\ (13) 
Jo 

This transformation has the advantage of being able to reproduce the two 
main components of the scattering amplitude valid respectively in the for- 
ward and the backward scattering directions. The comparison between the 
modified scattering intensity and the Mie theory, will be given in the next 
section. 



4 Numerical results 

In this section, we will compute the scattering intensity i{6) = for 
homogeneous, large spheres {ka > 20), from the original eikonal model (0) 
and the modified one (|13[). The model is in fact valid for ka > 10, and 
according to a numerical point of view, the new version is strictly equiva- 
lent to the original one when the expression (H) is used. The figures (2-6) 
show the results obtained when the refractive index has an imaginary part 
sufficiently large (according to the relation written above). We observe that 
the modified eikonal formula (dashed curves) is in much better agreement 
with the Mie theory (solid curves) for scattering angles greater than 40", as 
long as \Im{m)\ka > 1 and ka >> 1. Indeed, in this angular region, since 
the reflection becomes the main component of the scattering amplitude, the 
intensity predicted by the original eikonal model (dotted curves) lies system- 
atically above the Mie theory. We observe that no oscillations appear for 
angles more or less greater than 40° because the internal multiple reflections 
in the sphere are not included in the new formulation, leading to a smooth 
description of the scattered intensity clearly seen in Fig. 6. In Fig. 7, we 
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have considered the case of an elhpsoid with ka = 50 and kb = 75 in the z 
direction. Although no comparison can be made with respect to the Mie the- 
ory, there is again a difference between the modified eikonal model (dashed 
curve) and the original one (dotted curve) for scattering angles above 40°. A 
comparison with a sphere of the same ka and index (Fig. 1) shows that the 
scattered intensity for the ellipsoid is lower than the former in the backward 
direction, an effect partially due to more absorption. 

5 Conclusion 

We have derived a modified version of the eikonal model valid for large ab- 
sorbing spheres {\Im{m)\ka > 1, ka » 1), by imposing that the asymptotic 
expressions of the eikonal model amplitude in the forward and backward 
limits reproduce the diffractive scattering amplitude and the backscattered 
amplitude obtained from geometrical optics. We have shown that one has 
to replace the coefficient j{9) in the formula (^) by the Fresnel reflection 
coefficient r±{9). Numerical results show a great improvement of the eikonal 
formulation, especially for scattering angles greater than 40° and a much 
closer agreement with the Mie theory. 

T. Lemaire thanks the Brazilian agency CNPq for his financial support 
(project 300792/93-0). 
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Figure Captions 

Fig.l Geometry and notations. 

Fig. 2 Scattering intensity for a sphere of size parameter ka = 50 and re- 
fractive index m — 1.33 + i0.05. Mie theory (sohd), modified eikonal 
model (dashed) , original eikonal model (dotted) . 

Fig. 3 Scattering intensity for a sphere of size parameter ka — 20 and refrac- 
tive index m— 1.5 -|- zO.l. Mie theory (sohd), modified eikonal model 
(dashed) , original eikonal model (dotted) . 

Fig.4 Scattering intensity for a sphere of size parameter ka — 60 and refrac- 
tive index m = 2.5 + i0.05. Mie theory (solid), modified eikonal model 
(dashed) , original eikonal model (dotted) . 

Fig. 5 Scattering intensity for a sphere of size parameter ka — 20 and refrac- 
tive index m — 5.0 + iO.l. Mie theory (solid), modified eikonal model 
(dashed), original eikonal model (dotted). 

Fig. 6 Scattering intensity for a sphere of size parameter ka — 100 and re- 
fractive index m = 2.5 + zO.Ol. Mic theory (solid), modified eikonal 
model (dashed), original eikonal model (dotted). 

Fig. 7 Scattering intensity for an ellipsoid of size parameter ka = 50, kb = 
75 and refractive index m = 1.33 + i0.05. Modified eikonal model 
(dashed), original eikonal model (dotted). 
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